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Abstract 

This  paper  is  an  initial  attempt  to  harness  the  power  of  adaptive  importance  sampling  techniques  for  estimating  false  alarm  probabilities 
of  space-time  adaptive  radar  detectors.  Fast  simulation  using  importance  sampling  have  been  notably  successful  in  the  study  of  conventional 
constant  false  alarm  rate  radar  detectors,  and  in  several  other  applications.  The  principal  task  here  is  to  examine  the  viability  of  using 
importance  sampling  methods  for  STAP  detection.  Though  a  modest  beginning,  the  adaptive  matched  filter  detection  algorithm  is  analyzed 
successfully  using  fast  simulation.  Two  biasing  methods  are  proposed  in  this  paper.  One  of  them  known  as  the  g-method  is  shown  to  yield 
excellent  results.  The  important  problem  of  detector  threshold  determination  is  also  addressed,  with  matching  outcome.  The  work  reported 
here  serves  to  pave  the  way  to  development  of  more  advanced  estimation  techniques  that  can  facilitate  design  of  powerful  and  robust  detection 
algorithms  designed  to  counter  hostile  and  heterogenous  clutter  environments. 


I.  Introduction 

Estimation  of  false  alarm  probabilities  of  space-time  adaptive  processing  (STAP)  detection  algorithms  is  examined 
here  using  forced  Monte  Carlo  or  importance  sampling  simulation  (IS) .  STAP  algorithms  are  of  considerable  importance 
for  radar  detection  in  airborne  and  spaceborne  applications.  They  are  notoriously  intensive  from  a  computational  point 
of  view,  with  the  more  advanced  (and  robust)  ones  being  also  analytically  difficult  to  quantize,  [2].  Therefore  it  is 
appropriate  to  attempt  to  develop  fast  simulation  methods  that  could  be  used  in  their  analysis  and  design.  Related 
work  by  [3]  considers  a  Chernoff  bound  approach  for  importance  sampling  applied  to  a  specific  class  of  STAP  algorithms. 

In  this  paper  we  use  lessons  learnt  from  developing  IS  techniques  for  characterizing  conventional  constant  false  alarm 
rate  (CFAR)  detectors,  [4],  and  describe  an  experiment  in  applying  them  to  STAP  detection.  The  starting  (and  ending) 
point  of  this  unpretentious  effort  is  the  celebrated  adaptive  matched  filter  (AMF)  derived  in  [1]  and  which  represents 
the  arrayed  version  of  the  workhorse  cell  averaging  CFAR  detector  for  conventional  radar  signal  processing  algorithms. 
The  false  alarm  probability  (FAP)  performance  of  the  AMF  detector  is  known  in  integral  form  and  can  be  numerically 
computed  to  any  desired  accuracy.  Thus  it  forms  a  suitable  basis  for  validating  our  simulation  experiments.  Two  specific 
IS  methods  (described  in  the  sequel)  are  presented  and  the  better  (and  also  easier)  one  is  implemented.  On  a  general 
note,  IS  is  the  chief  simulation  methodology  for  rare  event  estimation.  It  is  an  enduring  method  that  has  distinguished 
itself  in  several  areas  of  science  and  engineering.  Briefly,  IS  works  by  biasing  original  probability  distributions  in 
ways  that  accelerate  the  occurrences  of  rare  events,  conducting  simulations  with  these  new  distributions,  and  then 
compensating  the  obtained  results  for  the  changes  made.  The  principal  consequence  of  this  procedure  is  that  unbiased 
probability  estimates  with  low  variances  are  obtained  quickly.  The  main  task  in  IS  therefore  is  determination  of  good 
simulation  distributions  for  an  application,  either  as  a  one-shot  feat  or  adaptively.  Simulations  performed  using  such 
distributions  can  provide  enormous  speed-ups  if  they  are  chosen  with  due  care  and  mathematical  precision.  Indeed,  if 
applied  successfully,  simulation  lengths  needed  to  estimate  very  low  probabilities  become  (only)  weakly  dependent  on 
the  actual  probabilities.  It  is  thus  possible  to  evaluate  any  probability  in  reasonable  amounts  of  simulation  time. 

During  the  conduct  of  simulations  reported  herein,  some  issues  concerning  the  adaptive  IS  algorithms  used  arise,  and 
these  are  discussed  briefly.  More  investigation  is  required  into  them.  However,  the  positive  outcome  of  the  methods  used 
is  that  excellent  match  with  numerical  results  is  obtained.  The  succeeding  sections  provide  a  short  statement  of  the  AMF 
algorithm,  how  IS  biasing  can  be  performed  to  hasten  false  alarm  events,  description  of  the  so  called  (/-method  which 
is  a  conditional  IS  technique  developed  originally  for  studying  sums  of  random  variables  ([5]),  the  fast  algorithms  used, 
how  inverse  IS  can  be  used  to  estimate  (and  therefore  design)  detector  thresholds,  simulation  results,  and  a  concluding 
discussion. 


II.  The  AMF  detector 

In  a  radar  system  consisting  of  a  linear  array  of  Ns  antenna  elements,  a  burst  of  Nt  pulses  is  transmitted  and  each 
element  receives  as  many  return  samples  in  any  one  range  gate.  The  NsNt  =  N  samples  are  complex  (because  of 
/  and  Q  channel  processing)  and  are  referred  to  as  the  primary  data.  They  may  contain  a  target  and  represent  the 
range  gate  to  be  tested.  The  samples  are  arranged  in  an  IV  x  1  column  vector  and  denoted  as  x.  The  target  return 
is  modelled  as  consisting  of  a  known  direction  vector  s  with  an  unknown  complex  amplitude  in  addition  to  clutter, 
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interference,  and  noise.  There  are  L  other  IV-length  complex  vectors,  called  the  secondary  data,  obtained  from  as  many 
nearby  range  gates  and  assumed  to  be  free  of  target  signal.  These  are  denoted  as  x(l),  l  =  1, . . .  ,L.  The  primary  and 
secondary  data  vectors  are  assumed  to  be  jointly  independent  and  complex  Gaussian,  sharing  the  N  x  N  covariance 
matrix  R  =  E{XX*},  where  the  superscript  f  denotes  complex  conjugate  transpose.  Under  these  assumptions  the  AMF 
detection  test,  as  obtained  in  [1],  is  given  by 


Is^R  1x|2 

— -  e  V 

stR_1s  H0 


(1) 


where 

R=iX>(l)x(l)t 

i=i 


is  the  estimated  covariance  matrix  of  x  based  on  the  secondary  data  (also  referred  to  as  sample  matrix),  and  77  is  a 
threshold  used  to  set  the  FAP  at  some  desired  level.  This  test  has  the  CFAR  property.  The  FAP  a  of  the  test  is  known 
to  be  given  by 


a  = 


L\ 

(L  —  N  +  1)!(./V  —  2)! 


rl  XL  AT+l^j  _  X^N  2 

Jo  (1  +  11X/L)L~N+1 


(2) 


which  can  be  used  to  numerically  determine  the  threshold  setting  for  a  desired  FAP.  As  shown  in  [1],  the  test  in  (1)  can 
be  rewritten  as 


stR-ixl2 


H 1 

Ho 

r/s^R  1s 

= 

r/s^R_1RR_1s 

,  L 

= 

r/s^R-1—  y^x(I)x(I)tR~1s 

^  1=1 

L 

= 

|^yRr1x(l)x(l)tR-1s 

1=1 

L 

= 

|f>tR"1x(l)|2 
^  1=1 

(3) 

of  the  usual  CA-CFAR  test.  The  LHS  is 

a  square  law  detector,  being 

the  output  of  a  matched  filter  (matched  to  the  direction  s  in  which  the  array  is  steered)  for  incoherent  detection  using 
the  so-called  sample  matrix  inversion  beamformer  weights  R_1s.  The  RHS  represents  a  cell  averaging  term.  Further 
details  on  these  issues  can  be  found  in  the  references  mentioned  above. 


III.  False  alarm  probability  estimation  using  IS 

Two  methods  to  quickly  estimate  FAPs  are  two-dimensional  (2-d)  biasing  and  the  conditional  (/-method  procedure, 
described  in  this  section. 

A.  2-d  biasing 

To  estimate  FAP  using  IS,  we  make  the  following  observations.  Suppose  each  complex  sample  of  a  secondary  vector 
is  scaled  by  a  real  number  01/2.  This  has  the  effect  of  scaling  the  covariance  matrix  estimate  R  by  6.  Therefore,  as  far 
as  the  covariance  estimate  is  concerned,  both  sides  of  the  test  in  (3)  remain  unaffected  by  the  scaling.  However,  each 
secondary  vector  being  scaled  by  d1^2  results  in  a  scaling  of  the  RHS  by  6.  Hence  choosing  0  less  than  unity  will  have 
the  effect  of  compressing  the  density  function  of  the  random  threshold  of  the  test.  Further,  a  scaling  of  each  complex 
component  of  the  primary  vector  by  a  real  a1/2  will  achieve  a  scaling  of  the  LHS  of  the  test  by  a.  Thus,  choosing  a 
larger  and  6  smaller  than  unity  will  achieve  an  increase  in  the  frequency  of  occurrence  of  a  false  alarm  event  during 
simulation.  The  IS  optimization  problem  will  be  a  two-parameter  one. 

The  (unbiased)  IS  estimator,  using  (1),  can  be  expressed  as 


1  ^  ^ 
a  =  —  ^l(|stR-1X|2  >  r/stR-1s)W(X,XL;a,6>); 


/*  (4) 
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where  the  notation  ~  /*  means  that  all  random  variables  are  drawn  from  biased  distributions,  and  Xl  =  (X(l), . . . ,  X(L))T 
with  K  denoting  length  of  the  IS  simulation.  In  setting  up  their  joint  densities,  we  use  the  fact  that  the  FAP  of  the 
AMF  has  the  CFAR  property  and  is  independent  of  the  true  covariance  matrix  R.  This  is  true  under  the  assumption 
of  Gaussian  distributions  for  the  data.  In  such  a  case,  the  simulation  of  the  AMF  test  can  be  carried  out  for  data 
possessing  a  diagonal  covariance  matrix  I,  denoting  the  N  x  N  identity  matrix.  Therefore,  primary  and  secondary  data 
can  be  generated  as  complex  vectors  with  independent  components.  The  unbiased  joint  densities  are 


/(*)  = 


rN 


and  f  (xl)  = 


-Xix(l)tx(l) 


71' 


LN 


so  that 


/(x,xL) = 


-x+x-Ei  x(l)tx(l) 


7r(L+l)N 

With  scaling  performed  as  described  above,  the  biased  joint  density  takes  the  form 

e-lxW-fXi  xdTxCi) 

/*(x,XL)  = 


7r(L+l)NaN^LN 


resulting  in  the  weighting  function 


W (X,  XL;  a,  9)  4 


/(x,xL) 


/*(x,XL) 

=  CaN0LNeA'aeB'e 


where 


A  =  x^x,  B  =  ^x(l)tx(l),  and  C  =  e 


-(A+B) 


(5) 


The  variance  of  the  IS  estimator  a  can  be  expressed  as 


vara  =  —  [/(iz)  —  a2] 
K 


(6) 


where  v  is  the  vector  biasing  parameter  (a,  6)T  £  [1,  oo)  x  (0, 1].  Denoting  by  A  the  false  alarm  event  in  (4),  the  quantity 
/  is  given  by 


I(v)  =  P*{l(M)W2(X,XL;zz)} 

=  £7{1M)W(X,Xl;i/)} 


(7) 


where  the  expectation  A1*  proceeds  over  biased  distributions.  Minimization  of  var  a  with  respect  to  the  biasing  parameters 
is  equivalent  to  minimization  of  I  and  is  described  in  the  Appendix.  Although  not  implemented  here,  this  description 
has  been  included  since  it  is  foreseen  that  such  a  method  could  be  useful  in  situations  wherein  the  g-method  might  be 
difficult  to  apply. 

B.  The  g-method  estimator 

This  method  exploits  knowledge  of  underlying  distributions  more  effectively,  yielding  a  more  powerful  estimator. 
Additional  advantages  are  that  only  a  scalar  parameter  optimization  problem  needs  to  be  tackled  and  the  inverse  IS 
problem  (for  threshold  optimization  or  selection)  can  be  easily  solved.  The  FAP  can  be  written  as 


=  P(|stR-1X|2  >  ?7stR_1s|Ho) 

=  ^{PflstR^XI2  >  ??st R-1s|XL)  Hq)} 

=  £{s(xl)} 


(8) 


Note  that  the  conditioning  in  the  second  step  above  is  equivalent  to  the  condition  that  a  covariance  matrix  estimate  is 
given.  We  proceed  to  estimate  a  using  the  form  in  the  third  step  above. 

With  the  condition  in  mind  it  is  easy  to  show,  assuming  that  X  is  rotationally  invariant  and  Gaussian,  that  the  random 
variable  slR_1X  =  w^X  is  distributed  as  CA/"i(0,  w^Rw)  with  independent  real  and  imaginary  components,  and  the 
weight  vector  w  =  R_1s.  The  random  variable  Y  =  |slR_1X|2  therefore  is  exponential  and  has  density  function 


/fe|XL,H0)  = 


y/w*  Rw 

wtRw 


y  >  o 
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Therefore 

S(XL)  =  P(Y  >  77StR~1s|XL.  H0) 

^-f)sfR_1s/w^Rw 

Note  that  if  R  =  R,  then  </(Xl)  =  e”^  and  this  is  the  FAP  of  the  AMF  when  the  covariance  matrix  is  known.  As 
discussed  before,  we  are  simulating  with  homogeneous  data  possessing  an  identity  covariance  matrix,  that  is,  with  R  =  I. 
The  (/-method  IS  estimator  then  takes  the  form 


where 


1 

K 

1 


K 


£(/(XL)W(XL;0) 


1 

K 

=  *£• 


—rj  D 


W(XL;0); 


(9) 


D  ~  |w|2 

stR-^s 

st  (R— 1)2s 

Choosing  the  (single)  biasing  parameter  9  <  1  thus  produces  a  decrease  in  D,  thereby  causing  a  higher  frequency  of 
occurrence  of  the  false  alarm  event  or,  more  appropriately  in  this  case,  a  larger  value  of  the  function.  Note  that  use 
of  the  (/-method  obviates  the  need  to  bias  primary  data  vectors.  Determination  of  a  good  value  of  9  proceeds  as  before. 
The  weighting  function  is  simply 

kF(xL;  0)  =  #LNe-(1-1/0)B  (10) 

which  can  be  deduced  from  (5)  by  setting  a  =  1.  The  scaling  9  is  optimized  by 


9m+ 1  —  9m 


which  is  just  a  one-dimensional  version  of  (13).  Estimates  of  the  /-function  and  its  derivatives  are  given  by 

k 

I< 

1 


^.92(Xl)W2(Xl;0);  ~  f* 

1 

K 

X  E-92(Xl)W(Xl;  W(XL -,9); 


?'(0)  =  77  E/  .?2(XL)W (Xl;  0)Wge(XL;  9);  ~  f* 


1 

K 


See  Appendix  A  for  definition  of  the  above  quantities. 

B.l  Simulation  gain 

A  useful  measure  of  the  effectiveness  of  any  IS  procedure  is  the  simulation  gain  T.  It  is  the  ratio  of  simulation  lengths 
required  by  conventional  MC  and  IS  estimators  to  achieve  the  same  error  variance.  Setting  the  variance  in  (6)  equal  to 
(a  —  a2) / k  (being  the  MC  variance)  where  k  denotes  the  length  required  by  the  MC  estimator,  yields  the  gain 

rv  —  rv ^ 

r  = 


/( u)  —  a 2 

While  the  simulation  gain  is  useful  in  learning  how  much  faster  than  MC  an  IS  technique  is  in  terms  of  simulation 
length,  it  also  serves  the  purpose  of  comparing  different  IS  estimators.  In  actual  simulations,  an  estimate  of  T  is  made 
by  using  the  estimates  for  a  and  I.  The  (/-method  estimator  has  simulation  gain  given  by 

=  a -a2 
9  Ig(6 )  -  a2 


SUBMISSION  FOR  RADAR  2005  CONFERENCE 


5 


Fig.  1.  Numerically  computed  FAP  of  the  AMF  detector. 


where  Ig  =  E*{(72(Xl)W2(Xl;  0)},  and  it  can  be  estimated  during  simulation.  It  always  has  a  smaller  variance 
and  consequently  larger  gain  than  the  IS  estimator  discussed  in  the  previous  section.  Indeed,  without  IS  (W  =  1), 
Ig  =  E{g2(X. l)}  <  E{g(XL)}  =  a.  That  Ig  <  I  with  IS  was  proved  in  [4]  for  conventional  CFAR  detectors.  The  proof 
in  the  case  of  the  detectors  considered  here  is  similar,  and  will  be  omitted. 


C.  Inverse  IS  and  threshold  determination 

The  inverse  problem,  namely  that  of  finding  by  fast  simulation  the  value  of  detector  threshold  77  satisfying  a  prescribed 
FAP,  is  readily  solved  using  the  ^-method  estimator.  This  is  done  by  minimizing  the  stochastic  objective  function 

J(v)  =  [ag(v)  ~  ao]2 


where  a0  is  a  desired  FAP.  Minimization  of  J  with  respect  to  77  is  carried  out  by  the  algorithm 

_  .  c  ao  —  (XgiVm)  _ .  0 

Vm+ 1  —  rfm  T  dg  . — .  ,  in  —  1,2,... 

a'giVm) 

where  Sg  is  a  step-size  parameter  and  the  derivative  estimator  in  the  denominator  is  given  by 


a'g(Vm)  =  -^J2De~r'Dw(X^e^  ~f* 


(11) 


with  the  prime  indicating  derivative  with  respect  to  ij.  Note  in  passing  that  this  derivative  estimator  actually  estimates 
(negative  of)  the  probability  density  function  of  the  AMF  statistic  on  the  left  hand  side  of  (1)  under  Hq. 

IV.  Numerical  results 

The  FAP  a  obtained  by  direct  numerical  integration  of  (2)  is  shown  in  Figure  1  and  is  used  for  comparing  IS  results, 
which  are  displayed  in  the  remaining  figures.  The  AMF  detector  consists  of  L  =  704  secondary  vectors  each  of  length 
N  =  352.  Shown  in  Figure  2  is  one  instance  of  adaptive  IS  estimation  of  FAP  for  a  (known)  threshold  of  77  =  56.50432. 
Figures  3  and  4  depict  results  from  implementing  the  inverse  IS  algorithms.  These  are  estimated  threshold  settings 
and  associated  FAPs  respectively.  Optimum  biasing  parameters  are  shown  in  Figure  5  and  simulation  gains  obtained  in 
Figure  6.  It  is  evident  that  match  with  the  results  in  Figure  1  is  excellent  and  this  has  been  numerically  confirmed. 

A.  Discussion 

The  IS  simulation  results  obtained  here  appear  deceptively  smooth  and  certainly  beg  an  obvious  question.  Indeed,  an 
artifice  has  been  employed  here  to  generate  them.  It  was  used  by  the  first  author  in  previous  work  on  capacity  estimation 
of  MIMO  channels  and  elsewhere,  and  found  to  be  extremely  useful.  In  conducting  rare-event  simulations  of  systems  that 
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Fig.  2.  Convergence  of  FAP  using  adaptive  IS  algorithms. 
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Fig.  3.  Threshold  optimization  for  AMF  detector  using  inverse  IS  algorithms. 


involve  signal  processing  operations  that  are  mathematically  complex,  there  are  two  principal  issues  that  contribute  to 
simulation  time.  These  have  to  be  dealt  with  effectively.  The  first  issue  concerns  the  rare  event  itself  whose  probability 
is  being  sought,  and  this  can  of  course  be  handled  by  suitable  IS  techniques.  The  second  concerns  the  computational 
intensity  that  accompanies  the  signal  processing.  The  two  are  not  unrelated.  In  the  case  of  STAP  detectors,  the  chief 
processing  burden  is  from  inversion  of  the  sample  matrix.  It  is  a  daunting  task  to  conduct  conventional  Monte  Carlo 
simulations  that  involve  several  millions  of  trials  to  estimate  low  FAPs,  with  as  many  matrix  inversions.  Assuming  that  a 
good  IS  scheme  can  reduce  the  number  of  trials  to,  say,  only  a  few  thousands  would  still  be  computationally  demanding 
(a  case  in  point  being  the  three  thousand  352  x  352  matrices  that  were  inverted  here).  This  is  the  point  at  which  IS 
departs  from  conventional  Monte  Carlo  in  a  subtle  but  important  manner.  It  is  almost  totally  useless  to  run  the  same 
random  variables  through  a  system  in  a  straight  MC  simulation.  With  IS  however,  much  can  be  learnt  by  repeatedly 
using  the  (same)  random  variables.  In  fact,  this  is  one  of  the  powerful  features  that  adaptive  IS  (and  inverse  IS)  can 
embed  into  complex  system  simulation.  But  how  does  an  IS  scheme  become  effective  in  the  first  place?  Assume  that  we 
have  a  biasing  scheme  that  promises  to  be  effective  once  the  parameters  of  the  biasing  distributions  have  been  optimized. 
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Fig.  4.  FAP  estimates  resulting  from  threshold  optimization  algorithms. 
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Fig.  5.  Optimum  biasing  (scaling)  parameters  for  IS. 

For  large  systems  (in  the  sense  of  number  of  random  inputs  involved),  running  truly  randomized  IS  algorithms  adaptively 
could  become  demanding  as  pointed  out  above.  If  system  performance  can  be  characterized  in  terms  of  certain  random 
‘metrics’  (we  use  the  word  with  a  slight  abuse  of  terminology),  then  these  metrics  can  be  pre-computed  for  a  given  set 
of  input  variables,  and  used  repeatedly  (which,  in  complex  systems  such  as  STAP  detectors  eases  the  computational 
burden)  in  adaptive  biasing  optimization  algorithms.  These  latter  algorithms  themselves  usually  require  no  more  than 
100  iterations  and  can  be  extremely  fast.  Resulting  IS  simulation  gains  can  be  simultaneously  estimated  and  these  tell 
us  whether  we  need  more  or  less  pre-generated  variables  to  achieve  certain  accuracies.  Adjusting  this  latter  number, 
biasing  and  system  parameter  optimization  (inverse  IS)  algorithms  can  be  run,  once.  Thus  there  is  an  initial  stage  of 
at  most  a  few  steps  during  which  gains  are  estimated  based  on  pre-computed  metrics  and  the  number  of  these  metrics 
is  adjusted.  All  this  is  not  as  complicated  as  it  appears.  Turning  attention  to  (9),  (10),  and  (11),  the  only  two  random 
quantities  (or  metrics)  that  are  needed  to  estimate  the  FAP  and  associated  detection  threshold  are  B  and  D.  This  is  for 
the  ^-method.  For  2-d  biasing  the  only  additional  quantity  required  is  the  norm  A  of  the  primary  vector,  defined  just 
after  (5)  and  this  adds  almost  nothing  to  the  computation.  It  turned  out  that  generating  K  =  3000  random  instances 
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Fig.  6.  Estimated  simulation  gains. 


of  B  and  D  was  certainly  an  overkill.  If  one  looks  at  Figure  6,  the  gain  provided  by  IS  for  estimating  a  =  10-b  is  about 
106.  From  usual  asymptotic  normality  arguments,  [4],  it  follows  that  about  100  optimally  biased  trials  are  sufficient 
to  guarantee  an  absolute  estimation  accuracy  not  exceeding  10%  with  95%  confidence.  For  2-d  biasing,  the  simulation 
gain  will  be  somewhat  lower  but  the  essential  advantages  of  the  method  above  remain.  That  is,  handling  a  few  hundred 
inversions  (once)  is  not  at  all  a  tall  order.  This  method  can  produce  such  an  avalanche  of  results  that  it  is  tempting  to 
think  of  it  (with  a  slight  stretch  of  imagination)  as  a  ‘turbo-IS’.  The  above  ideas  certainly  need  quantification  but  it  is 
beyond  the  remit  of  this  short  paper  to  delve  deeper. 

An  interesting  observation  comes  from  Figure  5,  which  shows  that  the  biasing  parameter  is  very  close  to  unity  and 
has  a  small  spread  despite  the  range  of  FAPs  considered.  The  implication  is  that  the  (one-sided)  density  of  the  metric  D 
has  small  variance,  presumably  owing  to  the  choice  of  L  and  N.  Smaller  values  of  these  constants  would  probably  lead 
to  larger  spread  of  biasing  parameter.  In  the  actual  adaptations,  a  small  value  of  the  step-size  parameter  was  used  to 
ensure  gradual  but  safe  convergence.  This  explains  their  apparently  slow  nature  as  seen  in  the  figures.  While  configuring 
results  for  a  suite  of  system  parameters,  only  the  first  adaptation  need  be  somewhat  long;  subsequent  adaptations  can 
be  much  shorter  as  they  pick  up  starting  or  initial  values  from  the  previous  one. 

V.  Conclusions 

In  this  short  and  admittedly  modest  paper,  we  have  made  a  small  inroad  into  the  use  of  adaptive  IS  algorithms  to 
characterize  a  STAP  detector.  The  AMF  was  used  as  example  and  results  have  been  very  good.  The  chief  reasons  for 
this  are  that  we  were  able  to  invoke  the  (/-method  and  inverse  IS,  find  a  suitable  biasing  strategy  that  could  be  easily 
optimized  adaptively,  and  find  a  way  around  the  difficult  task  of  inverting  large  matrices  several  times  (as  described 
above).  The  hope  is  that  applications  to  other  STAP  configurations,  such  as  normalized  AMF  and  those  that  handle 
non-homogenous  clutter,  will  also  meet  with  success.  But  this  remains  to  be  seen  as  we  are  certainly  not  in  position  to 
predict  what  subtleties  these  other  detection  algorithms  can  throw  up.  It  is  clear  that  IS  is  still  in  its  infancy,  especially 
insofar  as  use  for  characterizing  modern  detection  algorithms  is  concerned.  The  simulation  experiments  conducted  here 
have  raised  questions  (at  least  in  the  authors’  minds!)  that  need  to  be  answered  definitively. 

Appendices 

I.  Adaptive  algorithms  for  2-d  biasing 

The  /-function  is  estimated  as 

/(iz)  =  1^1(A)W2(X,Xl;Iz);  ~f*  (12) 

1 

and  its  minimization  can  be  carried  out  using  the  2-dimensional  adaptive  algorithm 

V m-\- 1  =  1 (l/m) 


(13) 
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Here,  S  is  a  step-size  parameter  used  to  control  convergence,  and  m  is  the  index  of  recursion.  This  is  a  stochastic  Newton 
recursion.  It  achieves  minimization  of  /  by  estimating  a  solution  of 

V/(iz)  =  (la  Ig)T  =  0 

where  Ia  =  dl(u)/da  and  Ig  =  dl(v)/d6.  The  estimate  of  the  Jacobian  J  (which  is  the  Hessian  matrix  of  I)  is  given  by 

j  _  /  laa  Ia6 

\Ia8  lee 


where  Ixy  =  dlx/dy.  It  is  straightforward  to  show  that  the  various  /-functions  defined  above  can  be  obtained  by  the 
notational  equations 


/,  =  E*{l(A)WWx} 
Ixx  =  Eic{l(A)WWxx} 
Ixy  =  E+{l(A)WWxy} 


with  various  derivatives  of  the  weighting  function  calculated  as 
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and  they  can  be  estimated  as  in  (12).  The  FAP  estimator  in  (4)  and  the  adaptive  biasing  algorithm  of  (13)  are 
implemented  simultaneously. 
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